Exact Scaling Relations In Relativistic Hydrodynamic Turbulence 
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We consider the steady state statistics of turbulence in general classes of dissipative hydrodynamic 
equations, where the fluctuations are sustained by a random source concentrated at large scales. It 
is well known that in some particular cases, such as non-relativistic incompressible turbulence, a 
Kolmogorov-type exact scaling relation for a correlation function holds. We show that all such scaling 
relations follow from a general relation on the current-density correlation function. The derivation 
does not require an energy cascade picture and suggests that this traditional interpretation of the 
Kolmogorov relation for incompressible turbulence may be misleading. Using this we derive exact 
scaling results for compressible turbulence in relativistic hydrodynamics, which reduce in the slow 
motion limit to the Kolmogorov relation. We discuss the experimental implications of the results. 
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Developed incompressible turbulence of fluids is one of 
the oldest problems of theoretical physics. Very little suc- 
cess has been achieved in the area due to the insurmount- 
able difficulties presented by the strong non-linearity. Es- 
sentially, the only non-trivial exact result on the statis- 
tics was derived in 1941 by Kolmogorov 0,0- He derived 
an exact relation stating that the third order correlation 
function of the flow velocity scales linearly with the dis- 
tance in a wide range of scales called the inertial inter- 
val. This relation is often interpreted as a statement on 
the energy cascade: the energy pumped into the system 
by a source at large scales cascades downscales due to 
flow instabilities. Eventually the energy is dissipated at 
a small scale by viscosity. In the wide intermediate range 
of scales, called inertial, viscous dissipation is negligible, 
the energy is conserved and is merely being passed on 
to smaller-scale fluctuations. Kolmogorov's relation ex- 
presses the constancy of the mean energy flux in the iner- 
tial range. The inertial range statistics is of much interest 
due to the conjecture that some properties of the statis- 
tics are universal in this range of scales. In particular, 
one expects a non-trivial scaling of the velocity correla- 
tion functions in the inertial range and the breakdown 
of self-similarity 0, Currently, the scaling problem 
seems inaccessible and the Kolmogorov relation is the 
only solid result. 

While it is natural to approach any problem of a turbu- 
lent type by first deriving a Kolmogorov-type relation for 
it, this was done actually only for quite a few problems. 
These are the passive scalar turbulence [1] , the magneto- 
hydrodynamic turbulence [f| Q and the Hall magnetohy- 
drodynamics |7j- An important example of a situation, 
where to the best of authors' knowledge no Kolmogorov- 
type relation has been derived is provided by the com- 
pressible turbulence. In this letter we describe a general 
relation for steady states of a turbulent type. This rela- 
tion reduces to known Kolmogorov-type relations in the 
situations listed above, while in other situations it gives 
new results. Our analysis indicates that the interpreta- 
tion of the Kolmogorov relation for the incompressible 
turbulence in terms of the energy cascade may be mis- 



leading, see also In fact, the relation can be found 
from the condition of stationarity of the pair-correlation 
function of velocity in the steady state, see e. g. @ and 
below, and no energy cascade picture is needed for the 
derivation. 

Using this we derive exact scaling results for compress- 
ible turbulence in relativistic hydrodynamics, which re- 
duces in the slow motion limit to the Kolmogorov rela- 
tion. We discuss the implications of the result, both in 
view of the possible applications to the description of the 
hydrodynamic behavior of the quark-gluon plasma and 
due to the existence of a dual gravitational description 
in terms of black hole geometry in asymptotically Anti- 
de-Sitter space. 

We consider a general class of classical field dynamics, 



d t q a + V-j a = f a , 



(1) 



where q a , a = 1, N are charges, j a are currents and f a 
are the external random source fields. These equations 
describe local conservation laws and provide a canoni- 
cal form for the effective hydrodynamic description. The 
latter studies the evolution of slow modes of the system 
provided naturally by the local conservation laws (since 
the zero wave-number component of charge fields is con- 
served, the low wave-number components evolve slowly 
by continuity [n|). The equations are closed via a con- 
stitutive relation that expresses currents in terms of the 
charges. In accord with the fact that hydrodynamics con- 
stitutes a low energy (wave-number) approximation, the 
relation has the form of a series in gradients, 



J, 



F?({P}) 



jb 



Gh b ({ P })V jP b 



(2) 



where dots stand for higher order terms involving either 
higher order derivatives or first order derivatives in power 
larger than one. Normally the zeroth order, reactive, 
term leads to a conservative dynamics while the first or- 
der term describes dissipation [lfj . For our purposes here 
the consideration of higher order terms in Eq. ^ is un- 
necessary and we will limit ourselves to the general class 
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of dynamics of the type, 



d t q a 
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The above truncation is usually sufficient for practical 
purposes. We will assume the standard mathematical 
formulation of the problem of turbulence where the forc- 
ing term f a is random and its statistics is stationary, 
spatially homogeneous and isotropic. This implies that 
the same properties of statistics hold for the steady state 
statistics of q a . The correlation length of the force will 
be denoted below by L. 

An important example of the above class is provided 
by the incompressible Navier-Stokes equations, 

V • v = 0, d t v + (v ■ V)w = -Vp + vV 2 v + f ,(4) 

where v is the flow velocity, p is the pressure, v is the 
kinematic viscosity and / is an external force field. It is 
for these equations that the Kolmogorov relation was de- 
rived originally jl[ . Consider the following version of the 
derivation, cf. [9j]. One considers the steady state con- 
dition dt{vi(0,t)vi(r,t)) — where the angular brackets 
stand for averaging over the force statistics. Using ho- 
mogeneity and isotropy it follows from Eq. (|U) 

= d t (vi(0,t)Vi(r,t)) = -2^(^(0, t)vi(r,t)vj(r,t)) 



+2( Vi (0, t)fi(r, *)) + 2vV 2 (v l {Q,t)v l {r,t)). 



(5) 



We now consider the limit of large correlation length L of 
the forcing. This allows us to consider the limit of r large 
but obeying r -C L. Due to r C L we have fi(r,t) « 
fi(0,t) and (Vi{0,t)fi(r,t)) « («i(0, <)/(((),«)) = e where 
e is some constant equal in this case to the energy pump- 
ing rate. At sufficiently large r the last term in the above 
equation is negligible, since it contains a higher number 
of spatial derivatives, and we obtain 

dj(vi(0,t)vi(r,t)vj(r,t)) = e, (6) 
(v i (0 ) t)v i (r,t)v j (r,t)) = er j /d, (7) 

where d is the spatial dimension and we used isotropy. 
The above relation on the triple correlation function is 
the Kolmogorov relation and using statistical symmetries 
it can be brought to the standard form 

([(v(r) - t7(0)) • r/r]) 3 ) = -I2er/d(d + 2) , (8) 

cf. 1. 

The general derivation of a Kolmogorov type relation 
for Eq. proceeds along the same lines. We consider 
the steady state condition dt(q a (Q,t)q a (r,t)) = (here 
there is no summation over a). Employing the dynamical 
equation ^ and using statistical symmetries, one finds 



= d t {q a (0,t)q a (r,t)) 
+2(q a (0,t)f a (r,t)) 
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Again we consider the limit of large correlation length 
L of the forcing. Keeping all parameters of the system 
fixed we study the limit where r in the above equa- 
tion is large but still much smaller than L. Because 
r «Lwe have f a (r, t) w f a (0,t) and {q a {0,t)f a (r,t)) w 
{q a (0, t)f a (0, t)) = e, where e is some constant. Taking 
the limit of large r we note that the last term becomes less 
important as containing larger number of spatial deriva- 
tives. Hence we obtain 



di{q a {0,t)F?{r,t)) = e. 
Assuming in addition isotropy one finds 



< 9 Q (0,i)^ a (r,£)) = ^i, 



(9) 



(10) 



where d is the space dimension. 

In the following we will consider the implications of 
this derivation to compressible turbulence in relativis- 
tic hydrodynamics. The implications to non-relativistic 
hydrodynamics are discussed in [8j . Consider first a rela- 
tivistic conformal hydrodynamics in (d + l)-dimensional 
space-time (d > 2). It is described by d + 1 fields: 
temperature T(x) and the (d + l)-velocity vector field 



M M (x),/i = 0, satisfying u^u^ 
energy tensor of the CFT obeys 



d v T^ v = 0, 







T. The stress- 



(11) 



and the equations of relativistic hydrodynamics are de- 
termined by the constitutive relation expressing T^ u in 
terms of the temperature and the four- velocity field. The 
constitutive relation has the form of a series in the small 
parameter (Knudsen number) 



Kn = l C or/L <C 1 , 



(12) 



where l cor is the correlation length of the fluid and L is 
the scale of variations of the macroscopic fields. Since the 
only dimensionfull parameter is the characteristic tem- 
perature of the fluid T, one has by dimensional analysis 



that l cor ~ -k. The constitutive relation reads 



oo 



T^(x) = > 2f», Tr - {Kn) 1 



1=0 



(13) 



where T^^x) is determined by the local values of ii M and 
T and their derivatives of a finite order. Keeping only the 
first term in the series gives ideal hydrodynamics, while 
dissipative hydrodynamics arises when one keeps the first 
two terms in the series. 

The ideal hydrodynamics approximation for does 
not contain the spatial derivatives of the fields. The I = 
term in (fT3|) gives the stress-energy tensor that reads (up 
to a multiplicative constant) 



T M „ = T d+1 [ Via , + (d+ l)u^u v ] 



(14) 



where rj^v = diag[- 
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The dissipative hydrodynamics is obtained by keeping 
the I = 1 term in the series in Eq. (fT3"|) . In the Lan- 
dau frame that fixes the ambiguity in the form of the 
stress-energy tensor under a field redefinition of the tem- 
perature and velocity, the stress-energy tensor reads (up 
to a multiplicative constant) 

= T d+1 [ Vliv + (d + l)u^u v ] - 2r)<x„ v , (15) 

where the shear tensor obeys a^u" = and is given 
by 

a ia> = 2 {^li u v + ®v u H + UvU P dpUn + Up,U P d p U v ) 

-^d a u a [r)p V + u^Uv] . (16) 

The dissipative hydrodynamics of a CFT is determined 
by only one kinetic coefficient - the shear viscosity r\. 

Consider next the addition of random force term at 
large scale L to (fTTj) . i.e. 

Q vT p, = p (17) 

Following the same procedure as outlined above, with 
the assumptions made about the steady state regime at 
length scales r much smaller than the force scale L much 
larger than the viscous scale I, one has the exact scaling 
relation 

£7* ' 

<T Oj -(0,t)T«(r,t)> = -^, (18) 

where (T O j(0,t)fj(0,t)) = e, and there is no summation 
over the index j. 

The exact result (fT8)l relies only on Eq. (fTT)) . Since 
the latter equation holds for general relativistic hydro- 
dynamics (also independently of the possible existence of 
additional hydrodynamic equations such as baryon num- 
ber conservation), then we conclude that Eq. (fT8|) is a 
general scaling relation for relativistic hydrodynamics. 
Analogous relation holds for (Tqo(0, t)Toi(r, t)) and cor- 
responds to the choice q — T o in Eq. (fTDj) . 

In the limit of non-relativistic macroscopic motions, 
the flow velocity v being much smaller than the speed 
of light c, the relativistic conformal hydrodynamics 
equations reduce to the non-relativistic incompressible 
Navier-Stokes equations [HI, [l3j]. Since the speed of 
sound in the relativistic conformal hydrodynamics is 
v s = -j=, the non-relativistic limit is the limit of small 
Mach number. In this limit, the exact scaling relation 
(fT8|) reduces to the Kolmogorov relation ([7|)- It is easy 
to see this explicitly, by expanding — (7,7^/0) and 
7 = [1 — v 2 /c 2 ]~ 1//2 in - and using the expansion of the 
temperature as 

T = T [l + P/c 2 + o(l/c 4 )] , (19) 

where P is the non-relativistic pressure. 

We can generalize further the above by considering 
charged hydrodynamics, where we have additional hydro- 
dynamic global symmetry currents conservation equation 

d»J» = . (20) 



By adding a random force / at large scale and the as- 
sumption of a steady state regime we have the exact 
scaling relation 

(J o (0,t)Ji(r,t)) = -j- , (21) 

where (Jo(0, t)f(0,t)) = e, and / is an external field. 

A particular interesting case is the hydrodynamics of 
anomalous relativistic gauge field theories, where the 
gauge invariant global symmetry current is conserved at 
the classical level, but is not conserved at the quantum 
level due to radiative corrections. The anomaly coeffi- 
cient C depends on the details of the microscopic theory 
and it is interesting to know how this coefficient is re- 
vealed in the hydrodynamic description [l4| [l9T |. In an 
external background electromagnetic field the anomalous 
equation reads 

dpT^ v = F» X J X , 9 M J" = CE^Bp , (22) 

where C is the anomaly coefficient and — 
F^u^B" = ^e^PuvFup. If we apply the same pro- 
cedure as above where now the external electromagnetic 
field is a random external source, we see that the anomaly 
coefficient appears in the correlation function (f2"Tj) . Note, 
however, that in order to actually extract the anomaly 
coefficient we will probably need a knowledge of an ad- 
ditional higher point correlation function in the steady 
state scaling regime. 

Consider next the experimental implications of the 
above exact scaling relations. The Reynolds number Re 
for the relativistic hydrodynamics can be estimated by 
the ratio of the ideal to viscous stress-energy tensors. 
This gives 

TL 

Re ~ — , (23) 

where T is the temperature, L is the characteristic length 
scale, and rj/s is the ratio of the shear viscosity 7/ and 
the entropy density s. A particulary interesting exper- 
imental setup is that of relativistic heavy-ion collisions 
such as the RHIC program. For gold collisions at RHIC, 
the characteristic scale L is the radius of a gold nucleus 
L ~ 6 Fermi, the temperature is the QCD scale T ~ 200 
MeV, and 2 ~ J- is a characteristic value of strongly 
coupled gauge theories. With these one gets R e ^> 1 and 
one may expect an experimental realization of the steady 
state relativistic turbulence [15( . It is therefore of much 
interest to see whether there exists in the RHIC data an 
experimental signature of the scaling relation (|18j) . 

Much like the numerical evaluation of anomalous ex- 
ponents in the scaling regime of non-relativistic incom- 
pressible fluids, one can attempt to study the existence 
of a scaling regime of relativistic hydrodynamics via nu- 
merical simulations. It would be interesting to use the 
existing relativistic hydrodynamic simulations for heavy 
ion collisions, where one can vary the system parameters 
such as the nucleus size L, in order to study the scaling 
relations at high Reynolds number. 
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A large class of compressible and incompressible hy- 
drodynamic systems have a geometrical realization as the 
dynamics of black hole event horizons in asymptotically 
Anti-de-Sitter space [ij], EH • In these descriptions, 
the forces arise from adding additional bulk fields, such as 
a dilaton, to Einstein gravity. It would be interesting to 
construct the geometrical analog of the above derivation 
of exact scaling relations. 
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